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, Abstract 

. A straightforward algebraic derivation is carried out of light quark 

CN I propagator and effective action in the instanton-anti-instanton molecule. 

' Exact expressions are obtained which contain contributions of all quark 

Q>^ . modes. Possible implications of the results include chiral phase tran- 

I sition and correlation functions of mesons and baryons. 

& 1 Introduction 

>■ . The theory of light quarks in the instanton vacuum has a long history - 

^ ! see e.g. review papers [1-3]. Recently a new impetus to the instanton liq- 

^ I uid modelQ has been provided by lattice measurements of quark and gluon 

■ - - ' correlators in coordinate space [4] showing the important role of instantons 

in forming hadronic wave functions and mass spectra. The standard theory 
of the instanton vacuum is based on the assumption of quark zero-modes 
dominance. This assumption was questioned in [5], and it was argued that 
non- zero-modes of each instanton play an important role in addition to zero- 
modes. To our knowledge it was the first systematic treatment of higher 
modes although their role in specific effects and within various approxima- 
tions were considered in some earlier papers [6-8]. 

Here we wish to present a complete treatment with all modes included 
of the so-called instanton-anti-instanton molecule, i.e. the // pair "bound" 

^Sometimes also called instanton gas model 
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through hght quarks. This object plays an important role in the instanton 
model of the QCD vacuum. According to the standard scenario the chiral 
symmetry restoration occurs due to the onset of the molecular component in 
the instanton-anti-instanton ensemble [9-15]. The analysis of chiral phase 
transition is not the purpose of the present paper. Our main objective is the 
detailed calculation of the molecule including quark propagator and effective 
action. Tooled with these results one can reexamine the chiral symmetry 
restoration in more details. Another even more immediate implication is 
the study of interactions in mesons and baryons induced by correlated // 
pairs [14]. The corresponding effective Lagrangian is derived from the quark 
propagator. As we shall see the quark propagator in the quantity in which 
the effects of non-zero-modes are mostly propounced. 

The systematic study of // molecule was performed in 1979 by C.Lee and 
W.A.Bardeen [16]. The propagator and the effective action were constructed 
making use of the power series in two independent parameters mp and p/d, 
where m is the quark dynamical mass, p is the size of (anti-) instanton, d- 
the distance between the two centers. In the propagator the interference 
terms between zero-modes and continuum were taken into account in the 
leading order with respect to the above parameters. Using this propagator 
the quark effective action was constructed. Surprisingly enough it contains 
only zero-modes overlap. In other words it was shown that the propagator 
is more affected by non-zero-modes than the effective action. The overlap of 
zero-modes in the effective action describes logarithmic attraction between 
instanton and anti-instanton. 

The revival of the interest to the role of non-zero-modes on one hand 
and to the // molecule on the other calls for the reexamination of the prob- 
lem. Below we develop an approach different from that of [16]. In its spirit 
our formalism is close to the multiple scattering theory [17,13], or scattering 
theory in a system with confinement [18]. Along these lines we have found 
the exact solution of the problem with all modes included. The final answers 
for the quark propagator and effective action are given in terms of matrix 
elements involving individual (anti-) instanton zero-modes and individual 
propagators with zero-modes deleted. Unfortunately the corresponding in- 
tegrals may not be evaluated analytically and we estimate them by the orders 
of magnitudes with respect to the parameters p/d and mp. In this way we 
present the non-zero-modes corrections to the effective action evaluated by 
Lee and Bardeen [16]. To get accurate results one has to perform numerical 
calculations which is beyond the scope of this paper. 
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The outline of this work is as foUows. In Section 2 wc introduce nota- 
tions and derive the "two-centers" equation for the quark propagator in // 
molecule. Section 3 is devoted to the evaluation of the quark propagator 
in // molecule. We present an exact expression for the propagator and an 
approximate form in which higher modes contribution is accounted for in 
the leading order. We argue that the interference of zero modes with con- 
tinuum is generally of the same order as zero modes contribution. Section 
4 deals with fermionic effective action. Again we find an exact answer with 
the contributions of all modes included. In Section 5 we discuss estimates 
for the propagator and for the effective action based on the existence of two 
small parameters - the ratio of the insatanton size to the size of the molecule 
and the quark current mass. Finally in Section 6 we discuss our results and 
suggest various ways of their implementations. 



2 Basic definitions and the general approach 

As usual in the instanton physics we work in Euclidean metric. Also in line 
with the common lore we use the superposition ansatz for the // system 

A,{x) = - Ri) = Kix) + Alix - d), (1) 

i=l,2 

where i = 1 corresponds to / placed at the origin = Rj = 0), i = 2 
corresponds to / placed at R2 = Rj = d. In singular gauge we can write 

^ -|- p^j 

where p and Jl/ are the size and color orientation of the instanton. Anti- 

instanton is assumed to have the same size p (this is not crucial), and its 
field follows from (2) though the substitution fj ^ rj and x ^ {x — d). Next 
we introduce the free propagator Sq, the single (anti-) instanton propagar Si 
with i = 1 for 7 and i — 2 ior I, and the total propagator S in the field (1): 

Sq^ — {—id — im), 

5-1 = {-id- gA^^ - im), 

S-^ = {-id - gA-im). (3) 
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Note that our free propagator 5*0 includes the current quark mass m which 
is not the case in Lee and Bardeen paper [16]. 

The instanton propagator has a formal spectral representation over nor- 
mal modes 

5, = V ' (4) 

with eigenfunctions > satisfying the equation 

{-td-gA')\ui>= \^\ui> . (5) 

Of special role are zero-modes with Aq = which have negative chirality 

-f5\u^ >^ -\u^ >, (6) 

where index n = is omitted. Anti-instanton zero-modes > have positive 
chirality. Zero-mode has the following explicit form 

{U^{x))ci = ^(^2a;2)l/2(^a;2 _|. p2)3/2 [( ^ )^]ij^jc'^&c = 



7r(2x2)l/2(a;2+p2)3/2 

where c and i are color and Lorentz indices, e is antisymmetric tensor, Tx = 
Tij.Xij.-i Tn — ("^jO) '^^ — ("^^ ~0 with T being the Pauli matrix. For anti- 
instanton X is replaced by {x — d), (1 — 75) ^ (1 + 75), r^; — > r^^. 

The non-zero-modes part of the instanton propagator (3-4) is known 
explicitly in the chiral limit m ^ [19] and given by [19,6,7] 



n>0 
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= Go{x -y){l + ^r,r;)(7r.7r,)-V2- (g) 
X y 

p^i'Kx'Ky)-^''^ f+r^-y 1+75 . T+_yf 1-75 

where Go{x—y) is the free propagator of massless fermion (not to be confused 
with 5*0 given by (3) in which the quark current mass is included) 



and 

tt.^I + pV'^'- (10) 

The propagator Sj{x,y) obeys the Green's function equation with m — 
and zero-mode deleted 

{-id-gA^)S'j^-a-\u^ ><u^\. (11) 

Mass corrections of the order 0(m) to Sj can be easily found from the spectral 
representation (4) rewritten in the form 

In the chiral limit the second term in the right-hand side yields Sj (see the 
first line in (8)). To get 0(m) correction from the third term we note that 
due to the orthogonality of eigenmodes 



lul >< ull ,v-^ \ui >< ul^ 



\"'n "-nl _ I "-nl -^2 _ ^^'^2 

n>0 n>0 ''^^ 

Thus the quark propagator in the single instanton field up to order 
terms has form [6] 

Si{x,y) = : hSj{x,y)+ 

—im 

+ imj dzS'j{x,z)S'j{z,y). (12) 

We note in passing that the instanton hquid model (see [2,3] and refer- 
ences therein) relies on the following approximation for the quark propagator 
in the single instanton field 

Sj^^-^^^^^ + {-td)-\ (13) 
—im 

The validity of this ansatz in instanton liquid model was discussed in [5]. 
The situation for the // molecule will be elucidated in what follows. 

For further purposes it will be convenient to recast (12) into the following 
form 

Si=^- —^ + So + Ai, (14) 

— 
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so that 

Ai = Ai = S'j-So + tm{S'if = {S'j - Go) + tm[{S'jf - Gl] = 

= S,+im[{S',f-Gl], (15) 
where Go is given by (9). For I similar notations are introduced with Aj = 

Now we derive the two-centers equation for the propagator S defined by 
(3). It follows straightforwardly from the superposition ansatz (1) and the 
definitions (3). Namely one can write 

^ S 1 s - 



_ \ Q^^ _ Q^^ \ 

Sq — {Si — So)Sq {S2 — So) So — {S2 — So)So {Si — So) 

(16) 

Generalization for more than two centers may be found in [13]. Next task is 
the actual calculation of (16), i.e. the invertion of the operator in the denom- 
inator. This will be done exactly without resorting to the approximations of 
the type 

Sc^So + J2{S^-So). 



3 Quark propagation function in the molecule 

Using (14) the propagator (16) may be rewritten as 

S = S2Z-'Si, (17) 

Z^So- Ai5o-'A2 - ^.^"^ o-^A^ - AiS^'^""' ^< '''' 



-im —im 
+ \u^>Y^<u\ (18) 

Performing calculations one may use any of the two equivalent expressions 
(16) for S. At intermediate stages of calcutions the 1 <-> 2 permutation 
symmetry is lost as in (18). In final answers this symmetry is restored. 
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To find Z ^ we will use the technique based on the relation 

A + \i>a<n\ i+<^1ilC> ^ ^^1^' 

with A being an operator and a - c- number. Applying it to (18) one gets 

1 



—im —im 

K = So- A,Sq'A2. (23) 

To find W^'^ one has to apply equation (20) repeatedly twice. After some 
straightforward but tedious work of collecting all terms in (21) we arrive at 

Z-^ = K-' + '-II^K-^lu' X u'IR^K-' + '-II^^Ik-'R.Iu^ X u^\K-^- 
T42 + ^2 ^-i|^i ^2|^-i _ ^2i + A,i ^_,^^^^, u'\R,K-\ (24) 
where 

f3 =<u^\So^A2K-^\u^ >, 
=< u^\So'\u^ >,V2i =< m'I^o'V 
Ai2 =< M^IS'o-'Asi^^'AiS'o- V >, (25) 

X = ,5o - Ai5'o-^A2, K^So-A^S^'A,, 

D^m'^ + {V12 + Ai2)(y2i + A21) - im{a + P) - ap. 

Ri\u^ >= AiSo^\u^ >, < u^\Ri =< ii^'|-So~^Ai, i, j = 1, 2. 

The final step is to insert (24) into (17) and use (14) for S^. This yields the 
following expression for the propagator 

S = ^{So + A2)K-\So + Ai) + l{So + A,)k-\So + A2)+ 



^vm^^g^ + A2)ir->i X u^\k-\So + A2)+ 
^12 + ^12 ^ A2)i^->^ X u'\K-\So + Ai)- 



^"^ {So + Ai)i^->2 >< u^\k-\So + A2). (26) 



^21 + A21 



D 

This expression is an exact one. It relies only on the standard superposition 
assumption (1) for the molecular field and on the formal representation (14) 
for the individual instanton propagator. All individual (anti-) instanton 
modes are included. No any small parameter like p/d or mp was used in 
the derivation. To get from (26) an insight into the physical picture one 
should cither perform numerical calculations or resort to approximations and 
estimates. This last strategy will be used in what follows. 
Keeping in (26) only linear terms in Aj one obtains 

d 

- ^(1 + R2)\u' >< u^\{l + R,) - ^(1 + R,)\u'' >< u'\(l + R2), (27) 
where 

d = m^ + Vi2V2i-im{a + 13). (28) 

This approximation coincides with the result of Lee and Bardeen [16]. In 
the chiral limit, more precisely when <^ V12V21, one gets even simpler 
expression 

Sc^So + Ai + A2- V2{\1 + R2)\u' >< m'|(1 + i?i)- 

-yi2'(l + ^i)k' ><^^'l(l + ^2). (29) 

It is tempting to simphfy (29) further and neglect Aj and Ri (see e.g.[3,14]). 

However such approximation is not justified since as it was already mentioned 
in [16] and will be discussed below in Section 5 all terms in (29) are in general 
of the same order of magnitude. 
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4 Fermionic effective action 



We define the quark effective action in // molecule as 

We leave aside at least temporarily the questions of normalization and 
regularization - see, e.g. [20-22]. From (16) we get 

detS~^ = detS^^ ■ rfetS'g"^ ■ detZ (31) 

where the quantitity Z is given by Eq. (18) of the previous section. To 
calculate detZ use can be made of an equation similar to (20) 

det{A + IC > a < 77I) = detA ■ (1+ < 77|-^|C >). (32) 

It is a simple but tedious task to apply (31) repeatedly to (18). In fact no 
new calculations has be done since the needed quantity essentially coincides 
with D defined by (24-25). Precisely 

detZ = detSo ■ det{l - So^AiSo^A2)^. (33) 
Combining (30), (31) and (33) one obtains 

--{r-J - (^S^'^S^)"'''' - W.''^.)0 (34, 

If we now perform standard regularization and normalization procedures [20- 
22] the final result for the effective action would be 

Te// = riff + riff + Indetregil - S^' A^S^' A-')+ 

+ Nfln[m^ + {Vu + Au){V2i + A21) - im{a + fi) - a(3]. (35) 

This result is an exact one and does not rely upon any approximation ex- 
cept for the superposition ansatz (1). The first two terms describe effective 
actions of individual (anti-) instantons. The third term corresponds to no- 
zero modes contribution. Approximated by ^ — tr(5'o"^Ai5'o""'^A2) it exhibits 
exponential damping as compared to the last one. This last term is of the 
main interest. It consists of the individual zero modes contribution and the 
interference between zero modes and continuum. The well known approxi- 
mation to this term is Nfln{m'^ -\- V12V21) [16]. In the next section we shall 
estimate corrections to this expression using the exact result (35). 
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5 Approximations and estimates 



In Section 3 the fermion propagator (26) has been determined and in the 
previous section the fermion induced effective action (35) has been found. 
However both expressions are not easy to handle due to high-dimensional 
integrals which cannot be evaluated analytically. Therefore in the present 
section we shall resort to estimates based on the existence of two independent 
parameters supposed to be small. These are p/d and mp, where p is the 
(anti-) instanton size, d- the // separation, m - the current mass of the 
quark. The widely accepted values of p and d are: p ~ 0.3 fm, d ~ 1/m 
[1-4]. The key quantities entering into most integrals (matrix elements) are 
Aj given by (14-15). According to (15) Aj is splittcd into mass independent 
part 6i{x,y) = Sl{x,y) — GQ{x,y) and a contributions proportional to m. 
Therefore corrections with respect to p/d and mp are independent. 

Let us start by providing arguments in favor of the statement made after 
equation (29) for the propagator, namely that zero modes interference with 
continuum is of the same order as zero modes product. 

Consider the situation when both x and y are close to the instanton, 
x y p. Then using (7) for zero modes one arrives at the following 
estimate for zero modes product 

< x\u'^ >< u^\y >- p-^d-^. (36) 

In the same region the interference term Ri\u^ >< may be estimated 
using (8) and (9). This yields 

< x\Ri\u^ >< u^\y >~ ( j dz6i{x, z)dzU^{z — d)) < u^\y >~ 

~ (p^d-V^)p"' ~ P''d-\ (37) 

In (37) we neglected mass corrections and substituted Ai by Si and Sq^ by 
dz, then performed integration over the shell of size p around z r-' d which 
resulted in p^ volume. Also unjustified (at least in certain regions) is the 
neglect of Aj as compared to 5*0 in (29). For example when a; ~ ?/ ~ p and 
(x — y) ~ p both 5*0 and Ai are of the order p~^. On the other hand when 
X p and y ^ d zero modes product dominates over interference. 

Prom the above example it is clear that to make estimates one should 
consider different regions of integration, typically the "near" region x ~ y ~ 
p, the "far" region x y d, and the "mixed" region x p,y d. 
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Looking for the estimate with repect to p/d one substitutes Sq^ either by 
1/d, or by 1/p. Attention should be paid to chirahty arguments, namely the 
matrix elements < v}\'~i^'~iv...'~ip\u^ > and < u^l'^^'ju-'-lplu^ > are respectively 
nonzero for odd and even number of 7 matrices. 

Consider now the essential part of the effective action (35) 

f e// = Nfln[m^ + {Vu + ^12) (V^2i + ^21) - im{a + /3) - ap] (38) 

The overlaps of zero modes Vij{i,j = 1,2) are of the order Vij ~ p'^d~^ [1- 
3,16]. The leading correction to (38) in p/d stems from Aij while that in mp 
arises from a and (3. Let us begin with the correction in p/d. For definiteness 
consider ^412. We are interested in the leading order correction and therefore 
can substitute K (see (25)) by -S'o. Since we are looking for the correction 
with respect to p/d we substitute Aj by Si (see (15)). Thus the quantity 
under consideration is 

=< u^\S^^52Sq^5iSo^\u'^ > . (39) 

By chirality arguments 6i{i = 1,2) should not contain contributions propor- 
tional to (1 + 75)72 (formally such contributions are found in {S'^ — Go), see 
(8)). Next we write 

^12 / ^^{/ dxu^^{x)dj2{x - z)}dz{ j dydi{z - y)dyV?{y - d)] = 

dz.Uz)d,Jy{z). (40) 

Consider the contribution from the region x<p,zc::^d±p,yc:^d±p. 
In this region u^^d^ ~ p~^. To estimate 62{x — z) consider equation (8) 
for ^2 with all arguments shifted by d since 5*2 corresponds to / placed 
at the distance d from the origin. It is easy to conclude that the leading 
contribution to 52{x — z) stems from the first term of the second line of (8) 
and S2{x — z) ^ Go{x — 2;) ~ d~^. Therefore Jx{z) ~ p'^p~^d~^ ~ pd~^. In 
Jy{z) one has dyu'^{y — d) p~^ . The estimate of 5i{z — y) deserves some 
attention. It also receives the main contribution from the first term of the 
second line of (18) but this time an additional small parameter arises: 

5,{z - y) ~ Go{z - y)[(7r.7r,)-V2(i + ^r.r^) - 1] ~ 

X y 

~ Go{z - y)pH-^ ^{z- y)-^pH-\ (41) 
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Therefore 

Jy{z) ~ J dy{z-y)-^d-\ 

where at the final step the integration over y will yield the factor. Now 
in (40) acting on {z — y)~^ in (41) will result in {z ~ y)~^ ~ since both 
z and y are ~ (dip). Substituting everything back into (40) we obtain 

A[, ~ p^d-^ (42) 

Since V12 itself is ~ p^d~^, the correction to it is of the order 0{{p/d)^). 

Now we turn to mp corrections. The leading one arises from a and (3. 
Consider 

a' =< u^\So^AiSo^\u^ > . (43) 

Chirality arguments are important here. If we substitute both 5*0"^ by 
Gq^ and Ai by 61 we get a zero. Therefore one should consider contributions 
arising from " unslashed" piece of any one of the three operators inside (43) 
and "slashed" pieces of the two others: 

a' ~ —im < u^ISiGq^Iu^ > —im < u^IGq^SiIu^ > + 

+ im< u^\Go\S[^ - GDGo^Iu^ im < u^\Go^51Gq^\u^ > . (44) 
The last matrix element may be represented as 

a' ~ m J dz{ J dxu^+{ X — d)dxSi{x — z)}x 

X {j dydi{z - y)dyU^{y -d)} = m j dzLx{z)Ly{z). (45) 

First consider the contribution from the "far" region x,y, z ^ d ± p. Then 
u^+{x - d)d^ ~ 5i{x -z)^ (see (41)), L^(^) - p^p-^d-^ - pd-^. 
The same estimate holds for Ly{z) and thus 

a' ~ mp\pd-^f ~ mp^d-^ (46) 

One may similarly verify the estimate (46) for the "mixed" region y ~ 
d ± p, z < p. Therefore the rri^ term in the effective action (38) acquires a 
correction of the order 0{p^/d^). 

Thus, we conclude that the leading-order corrections to the effective ac- 
tion (38) are given by 

Feff = Nfln{m^[l + Q{pVd^)] + l^i2^i + {V^2 + V^2i)0(pVt^')}, (47) 
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and we remind that Vij ~ p^/d^. We see that the corrections to V12V21 
in the effective action are of order p^/dP" ~ (1/3)^. A careful look into 
various normalization factors, powers of tt, etc. enhances this result by an 
order of magnitude. Here however a caution is needed since we are not 
in a position to account for possible numerical factors resulting from the 
integration procedure. 

Finally we wish to add one more remark on our results (38) and (47) for 
Fe//. The value of the overlap Vij being proportional to jS" depends also on 
the relative color orientation of / and /. Maximal attraction is reached under 
the condition that the relative color orientation vector = tr(r2jr^r2]^) 
is parallel to [14]. On the other hand the matrix element A^j has a 
different dependence on the color orientation. Therefore the most favorable 
color orientations for VijVji and VijAji entering into Fe// are different. The 
problem of the polarization of the // molecule with higher modes included 
will be discussed elsewhere. 

6 Summary 

In this parcr we have obtained exact sohitions for the quark propagator and 
the effective action in the instanton-anti-instanton mocecule. Up to now a 
consistent treatment of the non-zero-modes has been lacking. Possible appli- 
cations of the results might be twofold. The first one concerns the scenarious 
of the chiral symmetry breaking and restoration involving II molecules. De- 
tailed properties of the fermion determinant are relevant for this subject - see 
e.g. [11-15]. Another posible application involves quark interactions induced 
by // molecules and hadron correlators [14]. The corresponding studies were 
so far based on the quark propagator with non-zero- modes neglected. In 
such approximation the propagator exhibits the symmetry properties which 
are traceable in meson and baryon spectra. We have shown the relevance 
of non-zero-modes in the quark propagator. To which extent would higher 
modes affect the predicted properties of the spectrum? The resulrs of the 
present paper allow to study this problem. The work is in progress now. 
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